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Abstract 

New expressions for the multipole moments of an isolated post-Newtonian source, in the form of 
surface integrals in the outer near-zone, are derived. As an application we compute the "source" 
quadrupole moment of a Schwarzschild solution boosted to uniform velocity, at the third post- 
Newtonian (3PN) order. We show that the consideration of this boosted Schwarzschild solution 
(BSS) is enough to uniquely determine one of the ambiguity parameters in the recent computation 
of the gravitational wave generation by compact binaries at 3PN order: £ = —7/33. We argue that 
this value is the only one for which the Poincare invariance of the 3PN wave generation formalism 
is realized. As a check, we confirm the value of £ by a different method, based on the far-zone 
expansion of the BSS at fixed retarded time, and a calculation of the relevant non-linear multipole 
interactions in the external metric at the 3PN order. 
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I. INTRODUCTION 



The first aim of the present paper is to obtain, in the continuation of previous work 0, |^ , 
some new expressions for the multipole moments of an isolated gravitating source within the 
post-Newtonian (PN) approximation scheme of general relativity. The "source" multipole 
moments, as they were derived in Ref. (2], are given, for any multipolar order £, by certain 
functionals of the PN expansion of the pseudo-tensor of the matter and gravitational fields, 
formally valid up to any PN order. The moments have been obtained by solving the matching 
equation between the inner PN field and the outer multipolar expansion. Thev represent 
the generalization of the 1PN multipole moments obtained earlier in Refs. 0, Hi. 1 The 
source moments are the ones which parametrize the u multipolar-post-Minkowskian" (MPM) 
expansion of the external field of the source, as it was defined and investigated in Ref. 0]. 

In Ref. the PN moments have been given in the form of some volume integrals, whose 
support is non-compact because of the contribution of the gravitational field, and extending 
therefore over the full three-dimensional space. However, as we recall below, their integrand 
is made from the PN expansion of the stress-energy pseudo tensor t^ v , whose physical 
validity is restricted to the near-zone of the source. The formal extension of the integrals 
to r e |x| — > +oo was (uniquely) defined by introducing a suitable procedure of analytic 
continuation. 

In the present paper, we derive some equivalent expressions for the PN source moments, 
in the form of surface integrals, formally performed at r —>■ +oo. However let us emphasize 
that the physical meaning of the limit r — > +oo (which is taken after the PN limit c — > +oo), 
corresponds to considering what one can call the outer near-zone, namely the region which 
is at once far from the source, r ^> (source radius), and still within the near- zone, i.e. r <C 
(wavelength). For some physical problems, such alternative surface-integrals expressions of 
the PN moments can be quite useful, as the next result of our paper will illustrate. 

Indeed, our second aim is to make a contribution to the problem of theoretical templates of 
inspiralling compact binaries for gravitational- wave experiments such as LIGO and VIRGO. 
Current calculations of the gravitational waves are based on the expressions of the PN source 
multipole moments 0, 0| , as well as on high-order post-Minkowskian iteration of the external 
field 0,0,1. The mass-type quadrupole moment has been computed in the case of compact 
binary systems up to the 2PN order j^, 0] and more recently at 3PN order [ill . 12l |. The 



physical motivation for such high-order PN calculations of compact binary inspiral can be 
found in Refs. jH Q H H H H . 



It has been shown in that at 3PN order the radiation field of compact binaries, 
modelled as systems of point particles, contains three ambiguity parameters, £, k and (, due 
to an incompleteness of the Hadamard self-field regularization, used in the computation of 
the quadrupole moment of point particles at the 3PN order. In the present paper we shall 
show how one can determine the value of one of these parameters, (, without any use of a 
self-field regularization scheme. The idea is to consider the situation where one of the two 
masses is exactly zero. Such a limiting case corresponds to a single particle moving with some 
uniform velocity. As it turns out, the ambiguity parameter (, but only this one, survives 
in this limit. We can therefore compute it from the particular case of the 3PN quadrupole 
moment generated by a single object, specifically a spherically symmetric extended matter 



1 As usual the nPN approximation refers to all the terms up to the relative order ~ l/c 2n , where c is the 
speed of light. Powers of c and the gravitational constant G will generally be explicitly displayed here. 
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distribution, moving with uniform velocity in the preferred reference frame with respect 
to which the multipole moments are defined. The external gravitational field of such an 
object is evidently physically equivalent to a boosted Schwarzs child solution (BSS), i.e. a 
Schwarzschild solution viewed in a frame which is obtained from the usual "Schwarzschild 
rest frame" by a Lorentz boost. 

To compute the multipole moments of the BSS we propose and implement two methods. 
The first one is to insert the metric of the BSS into our new expressions for the PN mul- 
tipole moments in terms of surface integrals in the outer near-zone. The second method, 
substantially more involved, consists of computing the far-zone or "radiative" moments 2 of 
the BSS by expanding the metric in the far zone at fixed retarded time, and comparing the 
result with an analytic calculation of the various non-linear multipole interactions occurring 
in the far-zone moments up to the 3PN order. The source-type moment will thereby be 
determined in an indirect way. We find that the results of the two methods agree, and 
uniquely determine 

( = -—. (1.1) 
S 33 V ; 

We argue below that (jl.lj) is the unique value for which the 3PN wave generation formal- 
ism, as applied to binary systems in Ref. 0, 2|> incorporates the global Lorentz- Poincare 



invariance of general relativity. We have already reported elsewhere |20| that dimensional 
regularization leads to the same value for ( (as well as the determination of the values of 
the two other ambiguity parameters £ and k), and therefore has the important feature of 
"automatically" preserving the Poincare invariance of the formalism. The calculation of ( is 
similar to the one of its analogue in the 3PN equations of motion of point particle binaries, 
namely the so-called "kinetic" ambiguity constant uik which has also been fixed from 

the requirement of Poincare invariance, either by using an appropriate Lorentz-invariant ver- 
sion of Hadamard's regularization at the level of the equations of motion |23t l24j . or by direct 
imposition that the ADM Hamiltonian be compatible with the existence of phase-space gen- 
erators satisfying the Poincare algebra 25[. [The second and last ambiguity in the 3PN 



equations of motion is the "static" ambiguity constant u s [21, 22] equivalent to the am- 
biguity parameter A [i^, 24j[. Both ui s and A (as well as ui^ in fact) have been obtained 



by means of dimensional regularization |26l. |27|. The same results have also been achieved 
in Refs. jicl 129^ by expressing the equations of motion at 3PN order in terms of surface 
integrals surrounding the compact objects.] 

The paper is organized as follows. Section |Il] is devoted to a general investigation of the 
PN source multipole moments of extended matter sources. In Section Til Al we present some 
useful reminders of the formulation of the PN moments, and in Section Til Bl we obtain our 
new expressions of these moments in terms of surface integrals in the outer (or far) near-zone. 
Section UTTl is devoted to the investigation of the multipole moments of the BSS. In Section 
1111 Al we apply the new expressions of the PN source quadrupole moment to the case of the 
BSS and obtain the crucial coefficient which yields Eq. (|l.lj) by comparison with the BSS 
limit of the quadrupole moment of compact binaries. In Section 1111 Bl we expand the field 
of the BSS at retarded infinity and obtain the radiative-type quadrupole moment, which we 
then relate to the source-type moment by computing the non-linear multipole interactions 



2 Although the BSS does not radiate gravitational waves, it is sometimes convenient to employ a language 
for the different types of multipole moments which is similar to the one used in the more general case of 
non-stationary matter systems. 



3 



therein up to 3PN order (the technical details of the non-linear iteration are relegated to 
Appendix El). In this way we are able to confirm the value (jl.lj) . 



II. MULTIPOLE MOMENTS AS SURFACE INTEGRALS 



A. Reminders of the PN multipole moments 



As gravitational field variable we employ the standard "Gothic metric deviation" from flat 
space-time, that we shall subject all over this paper to the condition of harmonic coordinates, 
meaning that 

hT = y/^gv-rr, (2.1a) 
d v hT = 0, (2.1b) 

where g denotes the determinant and g^ u the inverse of the covariant metric: g = det(g pcr ) 
and g^gua — 8%, and where if" is the Minkowski metric written in Cartesian coordinates: 
j^nv _ diag(— 1, 1 ; 1 ; l). The Einstein field equations, relaxed by the harmonic-coordinates 
condition, read as 

□ w»> = T ^ ; (2.2) 

where □ = is the flat space-time d'Alembertian or wave operator, and t^ u denotes 
the total stress-energy pseudo tensor of the matter and gravitational fields in harmonic 
coordinates, given by 

r^ = \g\T^ + ^-^A^[h]. (2.3) 

Here T^ v is the matter stress-energy tensor, and A Miy is proportional to the stress-energy 
distribution of the gravitational field — K^ v is a functional of the field strength h or, more 
precisely, a function of h and its first and second space-time derivatives (~ dh and d 2 h), at 
least quadratic in h, dh or d 2 h. The pseudo tensor r^ v is conserved in a Minkowskian sense, 



d v T^ = , (2.4) 

which is equivalent to the equation of motion of the matter source: V„T M1/ = 0. 

In principle one can solve Eq. (|2.2|) by successive iterations in the order of non-linearity, 
i.e. by a formal expansion in powers of Newton's constant G. This would constitute what 
one calls the post- Minkowskian expansion (with explicit consideration of the matter source 
terms). However, such a straightforward post-Minkowskian expansion does not lead to eas- 
ily implementable iterations (see, e.g., [2O,0|). Another scheme, technically more useful, 
consists of splitting the problem of solving (|2.2|) into three sub-problems. First, one solves 
the vacuum equations, obtained by setting the matter stress-energy tensor to zero in 
Eq. (J2.3j) . by a particular post-Minkowskian expansion. Second, one solves the full inho- 
mogeneous Einstein equations (J2.2)) by a formal PN expansion, in powers of the velocity of 
light c (with G fixed). And third, one combines together these two expansions by means of 
an appropriate variant of the method of matched asymptotic expansions. 

The particular post-Minkowskian scheme we use to solve the vacuum Einstein equations 
is the so-called " multipolar-post-Minkowskian" (MPM) expansion of the metric exterior to 
the source 0. This expansion combines a non-linearity expansion in powers of G, with a 
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multipolar expansion of the successive non-linear iterations of h**", say hf^, i.e. essentially 
a decomposition of each function h^ u (t,r,8,ip) in tensorial spherical harmonics on the unit 
sphere parametrized by 9 and (p. It will be convenient, following 0, to denote the MPM ex- 
pansion of some quantity like h^, as M.{h lLV ), where the calligraphic letter Ai serves mainly 
the purpose of reminding the multipolar aspect of the expansion. [The post-Minkowskian 
aspect, though playing a crucial role in the definition and especially the construction of 
MPM metrics in |6J, is less important to keep in mind in the reasoning that we shall follow 
below.] The MPM expansion of the metric is thus written as 



M(h 



/1U\ 



+00 

E 

71=1 



(2.5) 



where one must always have in mind that each term hl^ u (t,r,9,ip) is decomposed in tenso- 
rial spherical harmonics (for instance STF tensorial harmonics like in 0). The successive 
post-Minkowskian coefficients are constructed iteratively from the linearized approxi- 
mation hi" by solving the harmonic-coordinates vacuum field equations. These consist of 
the harmonicity condition d u M.(h^ u ) = and 



UM{hT) =M{A 



)1V\ 



(2.6) 



where M.{K^ U ) means in which we substitute for h, dh, d 2 h their multipolar expansion: 
_M(A Ml/ ) = K^ u [M.{h)\. Again we stress that A1(some field) refers to a quantity which is 
supposed to be decomposed in a spherical-harmonics expansion (with coefficients being some 
functions of t and r). 

The algorithm to explicitly construct the successive post-Minkowskian coefficients h£ u , 
from the knowledge of the linearized approximation , was given in Ref. . As everything 
depends on hi", the parametrization of the linearized approximation will determine the full 
MPM metric ([2.5)1 . In our approach the source multipole moments are defined as the "seed 
moments" that one introduces at the start of the MPM scheme to parametrize the linearized 
approximation h^ v . 

With full generality one can write the linearized metric hf, which satisfies □ h^ v = 
together with d u h^ u = 0, in the form 



K = Cm + *Vi + 0V - fM 



(2.7) 



where ipi is a linearized gauge transformation vector (satisfying □ ipl = 0, so the harmonic 
gauge condition is preserved), and where h^ nl represents a useful form of the linearized mul- 
tipolar metric, called "canonical" and introduced in Ref. j32(|. The mass-type and current- 
type source multipole moments, Ii.it) and Jiit), respectively, are symmetric and trace-free 
(STF) tensors with respect to their I indices (£ is the multipolar order). They parametrize 
the canonical metric via the following definition, 
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3 In addition, the components of the gauge transformation vector 

Y L 



Our notation is standard 

(p^ can be parametrized by four other sequences of multipole moments, called Wl, Xl 

' i.e. 



and Zl, also being STF in their indices L, in the way specified by Eq. (4.13) in 
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The complete set of moments parametrizing the linearized approximation (|2.7|) with ()2.8|) 
and (|2.9j) . i.e. {Il, Jl, Wl, X l , Y l , Z l }, is collectively referred to as the source multipole 
moments. All these moments are defined in such a way that they admit a non-zero finite 
Newtonian limit, when c — > +oo. It is clear that the most important of these moments 
are the mass- type moment II and the current-type one Jl- Indeed the other moments, 
Wl, ■ ■ ■ , Zl, parametrize a gauge transformation and thus do not play any physical role at 
the linearized order (though they do play a role at the non-linear level). 

In Section IIIIBI we shall also recall the definition of two and only two sets of moments, 
named the "canonical" moments, denoted by {Ml, Sl}, which are physically equivalent to 
the complete set of six source moments {II, Jl, • • • , -Zzi, in the sense that they describe 
the same external gravitational field. However, following [2], we prefer to reserve the name 
of source moments to the set II, Jl, ■ • • , Zl because they are connected via some analytic 
closed form expressions to the stress-energy tensor of a PN source. 

The MPM algorithm computes sequentially any of the non-linear coefficients in Eq. ()2.5|) 
as follows 0, l33| . 4 Suppose that the first n — 1 coefficients hi, ■ ■ ■ , h n ~i, where hi is given 
by Eqs. ()2.7J) - (j2.9|) . have been constructed. We have then to solve, at the n-th order, the 
inhomogeneous wave equation, 

□ /C = Af, (2.10) 



whose source term is known from the previous iterations, 



i.e. 



n 



i- 



j h n -i\, 

and where as always we also have to satisfy the coordinate condition d u h^ u = 0. The solu- 
tion, satisfying a condition of "stationarity in the past" ensuring that the correct boundary 



L = ii ■ ■ ■ ii denotes a multi-index composed of £ multipolar indices i\, ■ ■ ■ , ii\ 8l = ■ ■ ■ di e means a 
product of £ partial derivatives 9j = d/dx l ; similarly xl = ■ ■ -Xi e is a product of £ spatial vectors 
Xi = x l ] symmetric-trace-free products are denoted with hats so that £l = STF(xl); sometimes we 
shall also use some brackets surrounding the STF indices: x/£,\ = Xl] the dots refer to the partial time 
derivation; e a h% is the Levi-Civita totally antisymmetric symbol (such that £123 = 1); index symmetrization 
means (ij) = 13 • r = |x| and u = t — r/c. 
4 In this paper we adopt a slightly modified version of the MPM algorithm, defined in Section 2 of Ref . 33] , 
which is more convenient in practical computations. 
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condition at Minkowskian past null-infinity are satisfied, reads 

*c =<+<". (2.H) 

The first term represents the standard retarded integral operator (denoted D^ 1 below) acting 
on the non-linear source, but augmented by a specific regularization scheme to deal with the 
divergencies of the retarded integral introduced by the fact that the multipolar expansion 
diverges at the origin of the spatial coordinates, when r = |x| — ► [see e.g. Eq. (|2.8|) ]. 
Posing 

(2.12) 

we do solve the required wave equation, i.e. \3u% u = A£ u , provided that the finite part (FP) 
takes the meaning specified below. The second term in Eq. ()2.11j) represents a particular 
homogeneous solution, i.e. Ov£ u = 0, defined in such a way that the harmonic gauge 
condition d u h^ = is satisfied at this order (see Refs. HQ f o r the details). 

To define the FP process we multiply the source term in Eq. (|2.12j) by a factor (r/ro) B , 
where B £ C and r denotes some arbitrary length scale, we compute the -B-dependent 
retarded integral in the domain of the complex B plane in which it converges, i. e. for which 
$l(B) is initially a large enough positive number, and we define it in a neighborhood of 
the value of interest B = by analytic continuation. The finite part in ()2.12|) means the 
coefficient clq of the zero-th power of B in the Laurent expansion ^ %> B p (where p £ Z) of 
the retarded integral when B — > 0. We emphasize that the divergencies cured by the FP 
regularization in ()2.12|) are ultraviolet (U.V.) type divergencies (r — * 0). As we shall see the 
same FP regularization will be used in the multipole moments to deal with their infra-red 
(I.R.) divergencies (when r — > +oo). 

In order to prevent any confusion, let us clarify the meaning of the various expansions 
that we shall use, and of the limits r — > and r — > +oo that will arise in the present 
paper. First, we recall that any MPM-expanded quantity is always thought of as written 
as a double expansion: one expansion in powers of G and one in spherical harmonics, say 
h n = YIi^lFl, where L denotes a multi-index which carries an irreducible representation 
of the rotation group (we do not write space-time indices). Then, each of the coefficients of 
this double expansion is given by some explicit function of t and r, say Fi(t, r, c), where we 
have also indicated a dependence on the velocity of light c (cf. for instance the simple case 
of the linearized approximation above). An important technical aspect of our formalism is 
that we shall consider these functions Fi(t, r, c) in the whole range of the radial coordinate 
r, even if this range does not correspond to a space-time region where the corresponding 
expansion is physically valid. Moreover, we shall sometimes consider instead of the original 
MPM coefficients F^(t,r, c) their PN expansion (or "near-zone expansion"), which means 
technically a formal expansion in powers of 1/c (with possibly some powers of In c) . We shall 
denote the PN expansion of any quantity by an overline. For instance, F L (t,r,c) denotes 
the expansion in powers of 1/c of Fi,(t, r, c), when keeping fixed the variables r and t. 

We deal sometimes [as in Eq. (J2.12j) above] with functions Fi,(t, r, c) in the region r — > 0, 
which is mathematically well defined (by real analytic continuation in r) but which physically 
corresponds to a region where the vacuum MPM metric should be replaced by a solution 
of the inhomogeneous Einstein equations, and therefore where the actual physical function 
would be a different function of r (and t) than F^(t, r, c). On the other hand, the limit r — ► 
+oo is acceptable both mathematically and physically for Fi(t, r, c), because it comes from 
a post-Minkowskian expansion which is valid all over the exterior of the source. However, we 
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shall also mathematically deal with the PN re-expansion Fi(t, r, c) of the function -Fl(£, r, c), 
and then formally consider the function Fi(t,r,c) in the limit r — > +00. Mathematically, 
the behavior of Fi(t, r, c) when r — > +00 is again well defined (by real analytic continuation 
in r) for each term Fi(t,r,c) in the MPM expansion. Though the limit r — > +00 seems 
physically incorrect for a PN expansion, it is here technically (or mathematically) well- 
defined. We stress that the PN limit c — ■> +00 is taken before considering the r — > +00 
behavior. If we remember that the PN expansion is physically valid only in the near-zone 
of the source, defined as r C cT, where T is a characteristic time of variation of the source, 
we see that the physical domain of validity of the r — ► +00 expansion of the PN expanded 
functions Fi{t,r, c) actually corresponds to the outer part of the near-zone, i.e. when r is 
much larger than the size, say a, of the source (multipole expansion), but still significantly 
smaller than a gravitational wave-length A ~ cT. We shall often refer to this domain as the 
far near-zone. 

As explained in Ref. in order to be able to define the behavior for r — > and r — > +00 
of all the coefficients F^t, r, c) appearing in the successive iterations of the MPM scheme one 
needs to make some formal technical assumptions: One must start with only a finite number 
of "seed" multipole moments, assume that they are infinitely differentiable functions of time, 
and that they tend to some constants in the infinite past. In our approach we assume that 
these requirements are initially satisfied, and we formally take, at the end of the calculation, 
a limit where these requirements are relaxed (so that we extend our results to an infinite 
number of moments, which are not necessarily past-stationary). 

The MPM expansion, sketched above, must be completed by an expansion scheme which 
covers the source. This is done by considering also a PN approximation for solving the inho- 
mogeneous Einstein equations. The PN scheme is a priori valid in the near- zone (r -C cT), 
while the MPM one is valid in the exterior of the source (r > a). The two domains of validity 
overlap in the exterior near- zone (a < r cT). One then imposes a "matching condition" 
which will enable us to determine the values of the multipole moments as functionals of the 
PN source. If we denote as above by h^ u the PN expansion the matching condition can be 
expressed as 

M{bT) =M(h^), (2.13) 

which says that the multipolar re-expansion of the PN metric agrees, in the sense of 
formal series, with the near-zone re-expansion (also denoted with an overbar) of the MPM 
metric M. {h^ u ). If we consider that a multipolar expansion is essentially an expansion in 
inverse powers of r when r gets far-away from the source, we can roughly summarize the 
matching equation (j2.13|) as saying that the far expansion (r — > +00, t = const) of the near- 
zone metric — L.H.S. of Eq. ([2.13)1 — coincides with the near expansion (r/cT — > 0, t = 
const) of the far (multipolar-expanded) metric — R.H.S. of (j2.13|) . The common general 
structure of both sides of (J27T3I) will be given in Eqs. (|2~2Tj) and (J2~22"j) . 

In Ref. 0] the PN source moments were obtained as functionals of the PN expansion 
of the pseudo-stress energy tensor defined by Eq. ([2.3)1 . namely . For the main source 
moments II and Jl (with any t > 2) we get 

h{t) = i FPJ d 3 x(^) B ^ dz{6 e (z)x L (T 00 + r*)(x,t + zr/c) 

4(2£ + 1) . . N „ dr° , 

d e+1 {z) x iL ——■ (x, t + z rj c) 



(£+1) (2£ + 3) ^ y ' cdt 



S 



2(2£+l) d 2 ^ 

5t +2 {z) Xij L -^7j( x >* + z r / c ) } > (2.14a) 



(£ + l)(£ + 2)(2£ + 5) t+AK ' vlJ c 2 dt 2 

J L {t) = - FPe ab{ie I d 3 x(-) B I dz{5 e (z)x L - 1)a T b0 (x,t + zr/c) 
c b=o j \r / J_ x 



2£+l df bc 1 

S e+1 (z)x L -i) ac -^(^,t + zr/c) > , (2.14b) 



(£ + 2)(2£ + 3) 



where we recall that xl means the STF product of £ spatial vectors, xl = STF(a;j 1 ■ • ■ x ie 
The other source moments, Wl, Xl, Yl, Zl, are given by Eqs. (5.17)-(5.20) in Ref. (2, 
We shall give below their new expressions in terms of surface integrals. See also f2j for a 
discussion of the conserved monopole and dipole moments (having £ = 0, 1). 

A basic feature of these expressions is that the integral formally extends over the whole 
support of the PN expansion of the stress-energy pseudo-tensor, i.e. from r = up to 
infinity. As already emphasized, the formal series is physically meaningful only within 
the near- zone. Therefore the integrals (|2.14|) physically refer to a result obtained from near- 
zone quantities only (in the formal limit where c — > +oo). However, it was found convenient 
in Ref. to mathematically extend the integrals up to r — > +oo. This was made possible 
by the use of the prefactor (r/r ) B , together with a process of analytic continuation in the 
complex B plane. This shows up in Eqs. (|2.14j) as the crucial Finite Part (FP) operation, 
when B — > 0, which technically allows one to uniquely define integrals which would otherwise 
be I.R. divergent, i.e. divergent at their upper boundary, |x| — > +oo. See Refs. [HIH for the 
proof and details. 

Since Eqs. (|2.14JI are valid only in the sense of PN expansions, the operational meaning of 
the auxiliary integrals in ([2.14)1 . with respect to the variable z, is actually that of an infinite 
PN series, given by 

pl +oo , o \ 2fc 

J ^dz5 e (z)r^,t + zr/c) = £ a k4 ^-J r^(x,t), (2.15a) 

(2£ + 1)!! 

where au e = - — — ^ — . (2.15b) 

(2A;)!!(2£ + 2fc + 1)!! v ; 



The expression of the function 5i(z), given in Appendix B of Ref. [3( as 

(2£+l)\\ 



8t(z) = 



|^p(l - z 2 Y such that ^ dz6 e (z) = 1 , (2.16) 



is useful when manipulating formal PN expansions such as (|2.15p . but will not be used 
explicitly in the present investigation. 



B. The multipole moments as surface integrals 

In this Section we derive an alternative form of the PN source moments ()2.14j) in terms of 
two-dimensional surface integrals. Such a possibility of expressing the moments, for general 
£ and at any PN order, as some surface integrals is quite useful for practical purposes, as 
we shall show below when considering the application to the BSS case. In keeping with the 
fact, just explained, that the "volume integrals" Eqs. (|2.14|) physically involve only near- 
zone quantities, the "surface integrals" into which we shall transform Eqs. (|2.14|) physically 
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refer to an operation which extracts some coefficients in the "far near-zone" expansion of 
the gravitational field, i. e. in the expansion in increasing powers of 1/r of the PN-expanded 
near-zone metric. Technically, as our starting point (j2.14j) is made of integrals extended 
up to r -> +00, our mathematical manipulations below will involve "surface terms" on 
arbitrary large spheres r = 1Z. All our manipulations will be mathematically well-defined 
because of the properties of complex analytic continuation in B. 

The basic idea is to go from the "source term", r^", to the corresponding "solution" h " , 
via integrating by parts the Laplace operator present in = □ h . From Eq. (j2.15j) 
we have 

/ d 3 ^r B x L f dz 6 t {z)*» (^T [ d^r B+2k x L UhT , (2.17) 

J J-i n k=0 \ c J J 

in which we insert □ = A — (^) 2 on the R.H.S., and operate the Laplacian by parts using 
A(r B+2k x L ) = (B + 2k) (B + 21 + 2k + l)r B+2k ~ 2 x L . In the process we can ignore the all- 
integrated surface terms because they are identically zero by complex analytic continuation, 
from the case where the real part of B is chosen to be a large enough negative number. 
[The complete justification of this is as follows. In the present formalism we are actually 
dealing with the MPM metric given by the non-linearity expansion (|2.5|) . and we are working 
with the calculation of some given finite post-Minkowskian approximation n. Then the PN 
expansion of the post-Minkowskian metric coefficient h£ u , namely h n , will typically diverge 
at infinity, but not more than a certain finite power of r, say N(n), depending on n and such 
that lim^oo N(n) = 00. Using h!^ = O (r N ( n A it is then clear that the all-integrated terms 
in question are zero when we choose initially $t(B) + 2k + £ + N(n) + 1 < 0, hence they are 
zero by analytic continuation in B.] Using the expression of the coefficients (j2.15b|) . we are 
next led to the alternative expression 

J d 3 xr B x L j dz 6 e (z)T fa ' = 

i ^_^B(S + 2£ + 4A ; + l)a M (-|) / d^r B +^ xjT" . (2.18) 

A remarkable feature of this result, which is the basis of our new expressions, is the 
presence of an explicit factor B in front of the integral. The factor means that the result 
depends only on the occurrence of poles, oc 1/B P , in the boundary of the integral at infinity: 
r — > +00 with t = const. At this stage it is useful to write down the expressions of the 
moments II and Jl we obtain by substituting (J2.18|) back into (J2.14|) . These are 

4(2£+l)(i? + 2f + 4A: + 3) / d \ 2k+1 f «B + 2 k -2~ 
(e + l)(2£ + 3) ak > e+1 {-cdt) J d * r XtLh 

2(2£+l)(B + 2£ + 4k + 5) ( d \ 2k+2 f j3 B+2k _ 2 , ^ 
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+00 ( , , \ 2k „ 

Jl = j^-g FP B ro B e ab{k (5 + 2£ + 4A; + 1) a kji I — J / d 3 ^r m -' z x L ^ )a h" 



{21 + 1)(B + 2t + 4k + 3) / d \ 2/c+1 /• r +2 ,_ 2a r 6 C 



(£ + 2)(2£ + 3) K '^ L \cdtJ 



a k>e+1 — / cfxr^ fc z x L „ l)ac h 



(2.19b) 



Let us proceed further. Thanks to the factor B we can replace the integration domain of 
Eq. (J2.18)) by some outer domain of the type r > 1Z, where TZ denotes some large arbitrary 
constant radius. The integral over the inner domain r < TZ is always zero in the limit B — > 
because the integrand is constructed from t^ v , and we are considering extended regular PN 
sources, without singularities. Now, in the outer (but still near- zone) domain we can replace 
the PN metric coefficients by the expansion in increasing powers of 1/r of the PN- 
expanded metric, which is identical to the multipolar expansion of the PN-expanded metric. 
This is precisely the quantity which was already introduced in Eq. (|2.13|) and denoted there 
M. (h^ u ) . Hence we have 



-1 

d 3 xr B x L j dz8i{z)T^ 

Vs(S + 2£ + 4ife + l)a M (-5-) / d^r B+2k -H L M(h^). (2.20) 



J -I 

^4 +00 / , \ 2fc 



16nG 

k=0 

We want now to make use of a more explicit form of the far near- zone expansion 

whose general structure is known. It consists of terms proportional to arbitrary powers of 

1/r, and multiplied by powers of the logarithm of r. More precisely, 

Al(^)(x,t) = ^^<>,t), (2.21) 

a, b 

where a can take any positive or negative integer values, and b can be any positive integer: 
a G Z, b G N. The coefficients ip^ u b depend on the unit direction n = x/r and on the coor- 
dinate time t (in the harmonic coordinate system). The structure (|2.21j) for the multipolar 
expansion of the near-zone (PN-expanded) metric is a consequence, via the matching equa- 
tion (J2.13|) . of the corresponding result concerning the structure of the near-zone expansion 
of the MPM metric ()2.5|) . which has been proved in Ref. by using the properties of the 
MPM algorithm for the iteration of the metric. [Again we stress the fact that the result has 
been proved at some arbitrary but finite post-Minkowskian order n; see Eq. (5.4) in Ref. 0. 
In the present paper we assume, following , that we are always entitled to sum up formally 
the post-Minkowskian series, and to view a result like ()2.21|) as true in the sense of formal 
power series.] 

As a side remark (which is not essential for the following), note that in Ref. Q the near- 
zone expansion of the MPM metric was viewed as an expansion in "ascending" powers of 
r/cT, namely it was written in the form 

W(x, t) = rP ( ln r ) q O n > *) > (2-22) 



11 



where p G Z and q G N. Evidently, the expansions ()2.22j) and (|2.21|) are equivalent. The 
only difference is that ()2.21|) was ordered in ascending powers of 1/r, i.e. in "descending" 
powers of r. Actually, both expansions are formal Laurent-type expansions, valid in some 
intermediate range of radii: a < r < cT. Their coefficients are related by a simple re- 
ordering of the exponents of r, 



0n^) = C(n,i). 



(2.23) 



In the following, we shall use the notation ^ 6 (n, t), corresponding to Eq. (j2.21|) . for the 
coefficients of these equivalent expansions. 

Inserting (J2.21j) into (|2.2(J|) . we are therefore led to the computation of the integral 



d xr 



B+2fc-2 



a,b Jn 



drr 



B+2k+l-a 



(In rf dnn L <p^ b (n,t), (2.24) 



where dQ is the solid angle element associated with the unit direction n (and hi = x^jr 1 ') 
The radial integral can be trivially integrated by analytic continuation in B, with result 



drr 



B+2k-\ 



- a (lnrf 



n 



-Y 

dB J 



B + 2k + i-a + l 



[2.25) 



Remember that we are ultimately interested only in the analytic continuation of such inte- 
grals down to B = 0. And as an integral such as ()2.25|) is multiplied by a coefficient which 
is proportional to B, we must control the poles of Eq. ()2.25|) at B = 0. Those poles are in 
general multiple because of the presence of powers of In r in the expansion, and the consec- 
utive multiple differentiation with respect to B in Eq. ()2.25|) . The poles at B = clearly 
come from a single value of a, namely a = 2k + £ + 1. For that value, the "multiplicity" 
of the pole takes the value 6 + 1. Here a useful simplification comes from the fact that the 
factor in front of the integrals in (|2.2()j) is of the form ~ B{B + K). In other words, this 
factor contains only the first and second powers of B. Therefore, only the simple and double 
poles 1/B and 1/B 2 in (J2.25|) can contribute to the final result. Hence, we conclude that it 
is enough to consider the values b = 0, 1 for the exponent b of In r in the expansion (|2.21|K 
To express the result in the most convenient manner let us introduce a special notation 
for some relevant combination of far-near- zone coefficients <y9^(n, t), which as we just said 
correspond exclusively to the values a = i + 2k + 1 and b = or 1. Namely, 



ttf^n,*) = a M -(2£+4k+l)<p 



7 2AH-£fl.o( 



n,t)+ l-(2£+4fc+l)lnr )<p% +l+11 (n,t) , (2.26) 



in which we have absorbed the numerical coefficient au/ defined by (j2.15bjl . [Notice that 
the coefficients (p^ u b depend a priori on the scale r .] With this notation we then obtain 

FP Br B (B + 2£ + Ak + l)a k/ [ d 3 x r B+2k - 2 % L M (h^) = 4tt ( ra L tf^) , (2.27) 

B ~° Jr>TZ 

where the brackets refer to the spherical or angular average (at coordinate time t), i.e. 

(fi L H?%) (t) = J ^nx*C(n,f)- (2.28) 
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As we can see, any reference to the intermediate scale 1Z has completely disappeared. The 
quantities ()2.28|) are integrals over a unit sphere, and can rightly be referred to as "surface 
integrals". These surface integrals will be the basic blocks entering our new expressions 
for the multipole moments. If we wish to physically think of them as integrals over some 
two-surface surrounding the source, we can roughly consider that this two-surface is located 
at a radius 71, with a < 71 < cT. Anyway, the important point is that, as we have just 
remarked, the surface integrals (|2.28|) . and therefore the multipole moments, are strictly 
independent of the choice of the intermediate scale TZ which entered our reasoning. 

Finally, we are in a position to write down the following final results for the source 
multipole moments (j2.19a|) and (j2.19b|) . expressed solely in terms of the surface integrals of 
the type (l2~2~51) . 

„2 +°° ( / J \ 2k 



2 too ( / i 
fc=0 k v 7 



4(2^+1) MfV v 



c 



(£ + l)(2£ + 3) \cdtj 

2(2£+l) ( d \ 

+ {£ + 1){£ + 2){2£ + 5) \Edl) 
.3 +°° r / a \ 2k 



n ijL *Lf2> > ( 2 - 29a ) 



Jl ~ Tn £ab ( 



4 G 

k=0 



1 — n v \ / 



2£ + 1 / d \ 



( d \ 2k+l ^ 

[^) <^-i>ac^ +1 )j. (2.29b) 



(£ + 2)(2£ + 3) \cdtj 

The other source moments, Wl, Xl, Yl and Zl, which parametrize the gauge vector given 
by Eq. ()2.9|) . admit similar expressions, which can be derived by the same method. For these 
we give only the results: 

cS V^/ 2£+l ( d V* Ji0 . 
WL = AG g{ (£ + l)(2£ + 3) UJ { ""^ -- } 



2/ • 1 ( d \ 2k 



2(£ + l)(£ + 2)(2£ + 5) 
-4 ±"r no _l 1 / ^ \ 2fe 



<^^ +2 >j, (2.30a) 



Xi ^ g U + l)(f + 2)(2. + 5) UJ <^*M+ 2 >} • ( 2 ' 30b ) 

J /J \2fe 

* = sbE - 

A:=0 k X 7 

3(2* + 1) ( d\ 2k+ \ h «, v 



(£+l)(2£ + 3) 



2(2^+1) / d V 



c 
4G 



(£+l)(£ + 2)(2£ + 5) 



/ d \ 2k+2 1 
(^J (^*M+2>j> (2-30c) 



f 2£ + 1 / <i \ 1 

^(■«E{ (< + 2)(2< + 3 ) <«l-D-«Sw>}- P-SOd) 
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III. APPLICATION TO A BOOSTED SCHWARZSCHILD SOLUTION 



A. Source quadrupole moment of the BSS at 3PN order 

As an application of our explicit surface- integral formulas (|2.29|) . we wish to compute the 
source-type multipole moments of a spherically symmetric extended body moving with uni- 
form velocity. Remember that our formalism assumes, in principle, that we are dealing with 
regular, weakly self- grav itating bodies. We expect, because of the nice "effacing properties" 
of Einstein's theory [34J, that our final physical results, especially when they are expressed 
as surface integrals like in (j2.29|) . can be applied to more general sources, such as neutron 
stars or black holes. Indeed, we are going to confirm this expectation in the simplest possible 
case, that of an isolated spherically symmetric body which is known, by Birkhoff 's theorem, 
to generate a universal exterior gravitational field, given by the Schwarzschild solution. We 
shall therefore apply our formulas to a boosted Schwarzschild solution (BSS). Actually, in 
order to justify our use of the BSS in standard harmonic coordinates, we must dispose of a 
small technicality. 

This technicality concerns the non-uniqueness of harmonic coordinates for the 
Schwarzschild solution, even under the assumption of stationarity (in the rest frame) 
and spherical symmetry. Indeed, under these assumptions, and starting from the usual 
Schwarzschild-Droste radial coordinate, say r$, the (rest frame) radial coordinate of the 
most general harmonic coordinate system, say r = k(rs), must satisfy the differential equa- 
tion (see, e.g., Weinberg jj^, page 181) 



d 
dr s 



2 2GM \ dk 



rs C* rs )dr s 



2k. (3.1) 



The "standard" solution of Eq. (JHHJ), which is considered in all textbooks such as |3J|, reads 
simply 

r = A ,standard (r5)=r5 _GM_ ( g_ 2) 

cr 

In the black hole case, the solution ()3.2|) is the only one which is regular on the horizon, 
i.e. when rs = 2GM/c 2 [as will be clear from Eqs. ()3.4|) below]. However, in the case of 
the external metric of an extended spherically symmetric body, regularity on the horizon is 
not a relevant issue. What is relevant is that the solution of the external problem (|3.1|) be 
smoothly matched to a regular solution of the corresponding internal problem. As usual, 
this matching determines a unique solution everywhere. In general, this unique, everywhere 
regular, solution will correspond, in the exterior of the body, to a particular case of the 
general, two-parameter solution of the second-order differential equation ()3.1|) . The latter is 
of the form 

/ r< \ r x 

[rs) , (3.3) 



fcg eneral (rs) = ^ ^ _ GM_^ + ^ fe 



where k 2 (r s ) denotes the (uniquely defined) "radially decaying solution" of Eq. ()3.1|) . and 
where C\ and c 2 are two integration constants. Indeed, when considering the flat-space limit 
of Eq. (J3.1|) . it is easily seen that there are two independent solutions which behave, when 
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rs — > +00, as r$ and r s 2 respectively. An explicit expression for the decaying solution is 



hirs) (2,2,4,^ , (3.1a) 



r s \ c^rs 



F(2,2,4,z) 



6 



2 + -(2-z)ln(l-z) 

z 



(3.4b) 



We can always normalize ci to the value C\ — 1. Then, with the above definitions, C2 has 
the dimension of a length cubed. By considering in more detail the matching of the general 
solution of the harmonically relaxed Einstein equations at the 2PN level (see, e.g., the book 
by Fock (36[, page 322), one easily finds that the second integration constant is of the order 
of c 2 ~ (GM/c 2 ) 2 a, where a denotes the radius of the extended body under consideration. 
It is also easily checked that the constant c 2 parametrizes, at the linearized order, a gauge 
vector ip\ of the form <p\ oc c 2 di(l/r), and can thus be referred to as a "gauge parameter". 
Comparing to the general multipole decomposition (|2.9j) . we see that this gauge parameter 
c 2 corresponds to the monopole (l = 0) in the gauge multipole sequence X L . 

Contrarily to the multipole moments of stationary sources, which are geometric invariants 
(and can be expressed as surface integrals on a sphere at spatial infinity), the "source 
multipole moments" defined in Ref. [2j (and re-expressed above as surface integrals over 
spheres in some intermediate region, a r cT) are probably not geometric invariants. 
They are useful intermediate constructs, which allow one to compute physically invariant 
information, but their definition is linked to the choice of harmonic coordinates covering the 
source. This means that the various "gauge multipoles" Wl, X L ,Y L , Z L will influence, at 
some non-linear order of the MPM iteration, the values of the two sequences of "physical 
multipoles" : II, Jl- Therefore, one should expect that, at some non-linear order, the physical 
multipoles II, Jl of a boosted general, harmonic-coordinate spherically symmetric metric 
will start to depend on the value of the gauge parameter c 2 . 

Here, we are only interested in computing the quadrupole moment 7y of a boosted general 
spherically symmetric metric. We shall see below that the index structure of iy will be 
provided by the STF tensor product of the boost velocity V 1 with itself, denoted V^V^ 
(assuming that the origin of the coordinates is at the initial position of the center of symmetry 
of the BSS). Therefore, any contribution to L t j coming from the gauge parameter c 2 must 
contain, at least, the factors c 2 and V^V^ , and also the total mass M. Taking into account 
the dimensionality of c 2 ~ (GM/c 2 ) 2 a, which is that of a length cubed, it is easily seen 
that there is no way to generate such a contribution to 7^. Therefore, we conclude that the 
source quadrupole moment of a boosted general, harmonic-coordinate spherically symmetric 
metric is strictly equal to the source quadrupole moment of a boosted standard harmonic- 
coordinate Schwarzschild solution, obtained by setting c 2 = (and c\ = 1) in ()3.3j) . i.e. by 
choosing the standard harmonic radial coordinate ()3.2|) . 

In the following, we shall therefore consider only such a boosted Schwarzschild solution 
(BSS) in standard form. We shall sometimes refer to the source of this solution as a black 
hole (though, strictly speaking, one should always have in mind some extended spherical 



Here 

F(a , ft7 ,, )=1+ ggi + a ' a+ . i y, +i ) g 

7 1! 7(7+1) 2! 

denotes Gauss' hypergeometric function. 
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star). For simplicity, we shall translate the origin of the coordinate system so that it is 
located at the initial position of the black hole at coordinate time t = 0. With this choice 
of origin of the coordinates all the current-type moments Jl of the BSS are zero. We shall 
concentrate our attention on the mass-type quadrupole moment iy, that we shall compute 
at the 3PN order. 

Le us denote by x M = (c t, x) the global reference frame, in which the black hole is moving, 
and by = (cT, X) the rest frame of the black hole — both x M and are assumed to be 
harmonic coordinates. Let x l (t) be the rectilinear and uniform trajectory of the (center of 
symmetry of the) BSS in the global coordinates x^, and V = (V 1 ) be the constant coordinate 
velocity of the BSS, 

V* = *?®-. (3.5) 
dt v ; 

The rest frame X^ is transformed from the global one x^ by the Lorentz boost (for simplicity 
we consider a pure Lorentz boost without rotation of the spatial coordinates) 

a? = A^(V) X\ (3.6) 

whose components are explicitly given by 



A° (V) 


= 7, 


(3.7a) 


A*o(V) 


V i 

= A°,(V) = 7-, 


(3.7b) 


A^.(V) 


J 7 + 1 c 2 ' 


(3.7c) 


with 7 




(3.7d) 



As explained above, we can assume that the metric of the BSS in the rest frame X M takes 
the standard harmonic-coordinate Schwarzschild expression, which we write in terms of the 
Gothic metric deviation H^, satisfying d u H^ u = 0. [We use the same conventions as in 
Eq. 1)2. ip . with upper case letters referring to quantities associated with the BSS rest frame.] 
Hence, 



, (1 + Sf) 3 

Y GM 



H m = 1 - \ W . ( 3 - 8a ) 

c 2 R 

H i0 = 0, (3.8b) 

(-12 Juf2 

H 13 = —^W NiN ^ ( 3 ' 8c ) 

where M is the total mass, R = |X| and N l = X 1 / R. A well-known feature of the 
Schwarzschild metric in harmonic coordinates is that the spatial Gothic metric if * J is made 
of a single quadratic-order term oc G 2 as shown in Eq. ()3.8cj) . The Gothic metric deviation 
transforms like a Lorentz tensor so the metric of the BSS in the global frame x^ reads as 

bT{x) = A^ p A u a H pa (K- l x) , (3.9) 

in which the rest-frame coordinates have been expressed by means of the global ones, i.e. 
X M (x) = (A _1 )^x y , where the inverse Lorentz transformation is given by (A _1 )^(V) 



Ifi 



A I/ /i (V) = A M I/ (— V). In our explicit calculations (done with the software Mathematica) we 
employ the BSS metric in exactly the form given by Eq. ()3.9|) . The only problem is to derive 
the explicit relations giving the rest-frame radial coordinate R and the unit direction N l as 
functions of their global-frame counterparts r and ri 1 , of the global coordinate time t, and 
of the boost velocity V 1 . For these relations we find 



R 



i + c 2 (y 



nil 

r 



1/2 



n 



7V* ® + ££(Vn 



1 + C 2(y - 1) (i) 2 - 2^(Vn) (1) + 7 2 ^ 



,1/2 



(3.10a) 
(3.10b) 



where (Wi) = V • n = V^vP is the usual Euclidean scalar product. The latter formulation 
of the BSS metric, Eqs. (|3.8|) - (|3.10|) . is well adapted to our calculations because we have to 
perform, when computing the source multipole moments, an integration over the coordinate 
three-dimensional spatial slice x £ M 3 , with coordinate time t = const, and this is easily 
done using the explicit relations ()3.10|) . 

However, let us notice that the BSS metric (in standard harmonic coordinates) is best 
formulated in a manifestly Lorentz covariant way as follows: 



( 



\ 



1 + 



GM 



GM 
c 2 r± 



v?u v 



G 2 M 2 



J 



-— n »n v , 
c A r z , 



(3.11) 



where is the time-like unit four-velocity of the center of symmetry of the BSS, where 
is the space-like unit vector pointing from the BSS to the field point along the direction 
orthogonal (in a Minkowskian sense) to the world line of the BSS, and where r± denotes the 
orthogonal distance to the world line (square root of the interval). The expression ()3.11|) is 
completely equivalent to (and more elegant than) the more "coordinate-rooted" formulation 
(I3.8|) - (|3.10|) . We shall employ it in a future investigation [37(. 

We compute the quadrupole moment Iy of the BSS, following the prescriptions defined 
by Eq. ()2.29af . To this end we first expand W when c — > +oo, taking into account all 
the c's present both in the expression of the rest frame metric ()3.8|) as well as those coming 
from the Lorentz transformation ()3.9|) - ()3.10|) . In this process the boost velocity V is to be 
considered as a constant, "spectator", vector. Note in passing that, in the present problem, 
the characteristic size a of the source at time t is given by the displacement from the origin, 
a ~ V t, where V = |V|, while the near-zone corresponds to r <C ct. Therefore, the far near- 
zone, where we read off the multipole moments as some combination of expansion coefficients 
ip^Kn, t), is the domain V t <^ r <^ ct. We have evidently to assume that V c for this 
region to exist. 

We then first get the near-zone (or PN) expansion of the BSS metric, h , by expanding in 
inverse powers of c up to 3PN order. Next we compute the multipolar (or far) re-expansion 
of each of the PN coefficients when r — > +oo with t = const. In this way we obtain what 
we have denoted by Ai[h^ Us j in Eq. (|2.21|) . In the BBS case it is evident that the far-zone 
expansion (|2.21j) involves simply some powers of 1/r, without any logarithm of r [indeed, 
see e.g. Eq. (I3~TU1) ]. 
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With AdCh ) in hand we have the coefficients of the various powers of 1/r, and we 
obtain thereby the various quantities defined by Eq. ()2.26j) . It is then a simple matter 
to compute all the required angular averages present in our formula ()2.29aj) and to obtain 
the following 3PN mass quadrupole moment of the BSS (the angular brackets surrounding 
indices referring to the STF projection), 



BSS 



mt 2 V {i V j) 



{ _9_ 83 V 4 

+ 14 ~& + 168 c 1 " 



507 V 6 
1232 ~& 



7 c 6 Vc 8 



(3.12) 



The first term represents the standard Newtonian expression, multiplied here by a bunch of 
relativistic corrections. (Recall that we have chosen the origin of the coordinate system at 
the initial location of the BSS at t — 0.) 

The last term in Eq. (J3.12j) . with coefficient C = 4/7, is the most interesting for our 
purpose. It is purely of 3PN order, and it contains the first occurence of the gravitational 
constant G, which therefore arises in the quadrupole of the BSS only at 3PN order. This 
term is interesting because it corresponds to one of the regularization "ambiguities", due 
to an incompleteness of Hadamard's self-field regularization, which recently appeared in the 
calculation of the mass-type quadrupole moment of inspiralling point particle binaries at the 
3PN order 11, The associated ambiguity parameter was called (, and was introduced 



as a factor of 3PN terms in the quadrupole moment having the form ~ m\ v[ l v^ or m\ v^v^ , 
where mi and m 2 are the two point masses, and v\ , v\ are their coordinate velocities. 6 The 
parameter ( represents the analo gue of the "kinetic ambiguity" parameter uj^ in the 3PN 
Hamiltonian of compact binaries 21, 2^| . It is now clear that ( can be determined from 



what we shall now call the BSS limit of a binary system, which consists of setting one of 
the masses of the binary to be exactly zero, say m 2 = 0. 

We have computed the BSS limit of the 3PN mass-type quadrupole moment of compact 
binaries computed for general binary's orbits in Refs. |lll . Il2l . We have also inserted for the 
position of the first body y\= v\t in order to conform with our choice for the origin of the 
coordinates. In this way we obtain 



L 



BSS limit 

ij 



j-2 (i j) 

m\ t v{ v\ 



9 vf 83 v? 

H H + 

14 c 2 168 c 4 



507 v\ 
1232 7 Q 



+ 



232 
"63~ 



G 2 m\ 



(3.13) 



The comparison of Eqs. ()3.13|) and ()3.12|) reveals a complete match between the two results 
if and only if we have the expected agreement between the masses, i.e. M = mi, and 
the velocities, v \ = V % (since the velocity of the body remaining after taking the BSS limit 
should exactly take the boost velocity), and the ambiguity constant ( takes the unique value 



c = - 



33 



(3.14) 



6 See Section X of Ref. 0] for the discussion on why and how to introduce the ambiguity parameters in 
the 3PN quadrupole of point particle binaries. 
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Our conclusion, therefore, is that the ambiguity parameter ( is uniquely determined by 
the BSS limit. Because of the close relation between the BSS limit with Lorentz boosts, 
it is clear that ( is linked to the Lorentz-Poincare invariance of the multipole moment 
formalism of Ref. 0] as applied to compact binary systems in 11, l3|. This link strongly 
suggests that the specific value (|3.14j) represents the only one for which the expression of 
the 3PN quadrupole moment is compatible with the Poincare symmetry. In other words the 
present calculation indicates that the Poincare invariance should correctly be incorporated 
into the laws of transformation of the source-type multipole moments for general extended 
PN sources as given by Eqs. (j2.14|) or ()2.29j) . though we have not verified this directly. 

Note that the situation regarding ( is the same as for the kinetic ambiguity parameter Uk 
in the 3PN equations of motion, whose value has also been uniquely fixed by imposing the 
Lorentz-Poincare invariance of the formalism I23l. l24l. l2q . (The other ambiguities £ and k 
in the binary's 3PN mass quadrupole moment |ll. Il2l| parametrize some Galilean invariant 
terms which cannot be derived from a Lorentz transformation. This can also be seen from 
the fact that £ and k are in factor of some "interacting" terms, depending on both masses 
mi and m 2 , and are thus out of the scope of the BSS limit.) 

Let us finally emphasize that we have obtained Eq. ()3.14|) here without using any regu- 
larization scheme for curing the divergencies associated with the self field of point particles. 
However, it has been shown in recent related works (2(1 13?| that this value is precisely the 
one derived in the problem of point particles binaries at 3PN order by means of the di- 
mensional self-field regularization, instead of Hadamard's regularization. (This shows that 
dimensional regularization is able to correctly keep track of the global Poincare invariance 
of the general relativistic description of isolated systems.) 



B. Retarded far-zone field of the BSS at 3PN order 

In this Section we present, as a further confirmation of our result, an alternative derivation 
of the crucial coefficient C = 4/7 in Eq. (|3.12j) . and hence of the kinetic-type ambiguity 
parameter ( = —7/33. This new derivation will be based on the expansion of the BSS 
metric (|3.8|) - (|3.9|) at Minkowskian future null infinity, r — > +oo with u = t — r/c = const. It 
relies on the complete identification of the metric ()3.8j) - (j3.9p . considered in the limit where 
r — > +oo with u = const, with the multipolar-post-Minkowskian (MPM) external metric of 
an isolated system, as given by the quantity M.{h^ u ) reviewed in Section Til Al above. This 
identification is justified by the fact that the general harmonic-coordinate MPM exterior 
metric defined by Eq. ()2.5|) smoothly matches, via the matching equation ()2.13|) . the inner 
PN solution of the inhomogeneous field equations inside the matter source (see Ref. fH). 
Indeed, so does the BSS metric (|3.8|) - (|3.9j) . as we have seen during the discussion on the 
general solution of the harmonic-coordinates condition (|3.1|) for the Schwarzschild metric, 
and our proof that the parameter c% does not contribute to the quadrupole moment. 

In this Section we thus have to consider the BSS metric as a functional of the six sets of 
multipole moments {II, Jl, Wl, Xl, Yl, Zl} parametrizing the linearized approximation 
to the MPM metric, namely h^ v defined by Eqs. (J2.7|) - ()2.9p . and then take into account 
the subsequent non-linear iterations, hijf , , • • • , whose computation is defined by the 
MPM algorithm of Ref. Q, as reviewed in Eqs. (|2.11|) - (J2.12|) above. Actually we shall work 
mostly with the so-called " canonical" MPM metric, h^ n , a functional of two and only two 
types of multipole moments called "canonical", mass- type moment and current-type Sl, 
instead of the six source moments II, Jl, • • • , Z L . Working with the metric h^ n simplifies 
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drastically the computation of the non-linear interactions. We shall justify below why, for 
our purpose, we can indeed use just this special case of a particular harmonic-gauge-fixed 
external metric, rather than the more general harmonic-coordinate external metric. Those 
metrics are geometrically equivalent, but we shall have to check (as in the case of the C2 
parameter in Section 1111 A|) that gauge effects do not modify the result we are interested in. 

In the present Section we shall have to compute the quadratic and cubic metric correc- 
tions, h^ n 2 and h^ n3 in the notation of Eq. (|2.5|) . for some specific multipole interactions. 
The canonical MPM metric is defined by 
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where the linearized approximation is given by the same formulas as Eqs. ()2.8|) but with the 
canonical moments Ml and Sl in place of the source moments II and Jl, and with all the 
gauge multipoles {Wl, X l , Y l , Z l } set to zero. In other words, 
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Then the non-linear metrics h^ nn [ML, Sl], for n > 2, are obtained from ()3.16j) by means of 
the same algorithm as before, explained in Eqs. (j2.11j) - (J2.12j) . in the case where the gauge 
vector 99^ = 0, together with the replacement {II, Jl} {Ml, Sl}- 

We have already alluded to the important point, proved in Ref. 0, that although the 
canonical metric ()3.15|) - ()3.16|) is simpler than our previous construction (|2.5|) - (j2.9j) . it is 
physically or geometrically equivalent to it, i.e. it describes the same physical matter system, 
provided that Ml and Sl are related to II, Jl, • • • , Zl by some specific relations of the 
type 



(3.17a) 
(3.17b) 

where Tl and Gl denote two non-linear functionals (at least quadratic in the moments) of 



M L = + J, W, X, Y, Z] , 

S L = Jl + Gl [I, J, W, X, Y, Z] , 



the original set of source moments II, Jl, 



Zl- However, if the use of only two sets of 



moments, Ml and Sl, is very useful when computing the non-linear multipole interactions, 
it remains that these moments have still to be related to the more "fundamental" source 
moments by Eqs. ([3. 17)1 . Indeed we know the analytic closed- form expressions of the source 
moments II, Jl, • • • , Zl (see Section Til A|) but similar formulas for Ml, Sl, valid to all PN 
orders, are not known to exist. The equations ()3.17|) need to be investigated anew for each 
specific cases. Fortunately, M L and Sl are "almost" equal to their counterparts II and Jl- 
Indeed we know that in the case of the mass-type moments for instance, Eq. (J3.17a|) when 
further PN expanded reads as 



M l = Il + ^5Il + o(\ 



(3.11 
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where 5Il denotes some correction term arising at order 2.5PN only (if necessary this term 
is given by Eq. (4.24) in Ref. |38j). The remainder in Eq. (j3.18|) is of order 3.5PN. The 
equation (|3.18|) shows that the 3PN term in the quadrupole moment of the BSS we are 
looking for — last term in Eq. (J3.12)) — will be the same for ly as for My. In addition to 
the relations (|3.17|) - (|3.18|) we must also take into account the possible effect of the coordi- 
nate transformation between the canonical metric (J3.15|) and the metric (J2.5)) . since as we 
mentioned, it is the latter which must be identified with the retarded far-zone expansion of 
the BSS metric (EPl) . 

Let us consider the general structure of the mass- type moments Ml (or II) in the case 
of the BSS. In the present problem there is only one vector which can be used to build the 
moment: Namely the boost velocity V" 1 , so the index structure of Ml must necessarily be 
made of the STF product V^) = ■ • • V it y [Indeed we recall that we choose the origin 
of the coordinate system to lie on the trajectory of the BSS, so we do not have at our 
disposal the vectorial separation between the origin and the BSS world line. It is clear 
that such a restriction is not physically crucial, however it simplifies the presentation and 
several arguments in this Section very much.] In addition, we readily see on dimensional 
grounds that V^l) must be multiplied by either u e , where u = t — r/c, or by the product 
(GM/c 3 ) or by (G M / c 3 ) 2 u e ~ 2 , and so on, and that each of the latter terms can be 
multiplied by some relativistic corrections of the type (V 2 /c 2 ) n up to any PN order. Here 
M denotes the BSS constant mass monopole, I = 0, or ADM mass. The general structure 
of the mass-type moment Ml (and also of II as well) at the 3PN order therefore reads 



where the notation & means that we have to add a term having the structure that is indicated 
next. In the quadrupole case, £ = 2, we recognize in the last (explicit) term of Eq. (J3.19JI 
the interesting form of the contribution to the ambiguity ( in the BSS limit. 

The mass multipole moment of the BSS varies with time typically like some u . This fact 
seems to be incompatible with the construction of MPM metrics in Ref. since it was 
assumed there, in order to implement this construction, that the matter source is stationary 
before some fixed finite instant — T in the remote past (and also that the coordinates are 
mass-centered in the sense that the dipole moment Mj is always zero). These assumptions, 
made in Ref. for purely technical reasons, imply in principle that all the multipole mo- 
ments Ml are constant before the date — T (i.e. for u < — T), and thus cannot a priori 
be applied to the physical situation of the BSS. Nevertheless, we shall admit in the present 
paper that we are allowed to use the construction and the results of |fj even in the case of 
the BSS. Indeed there has been several indications in our previous works using the MPM 
formalism, notably when the formalism was used for the computation of gravitational wave 
tails and "tails-of-tails" 0, H[, that the MPM expansion can in fact be applied to more 
general sources which have always been non-stationary, for instance an inspiralling compact 
binary formed by capture of two particles initially moving on some hyperbolic-like orbits. 
On the other hand, as we shall see our application of the MPM formalism to the external 
field of the BSS will yield some consistent result, which is independent of any initial instant 
— T and is in agreement with the result of Section fill Al This justifies a posteriori (to some 




(3.19) 
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extent) our expectation that the MPM formalism is still valid in the case of the BSS. Fur- 
thermore, we shall give in Appendix El an explicit proof that the integration formulas of the 
MPM formalism admit a well-defined limit when the multipole moments are continuously 
deformed into those of the BSS. 

We are interested in the quadrupolar contribution in the full non-linearity expansion 
(|3.15|) . as seen from retarded infinity, r — > +00 with u = const (this limit will be referred 
below to as X + ). We base our investigation on the time-time component (00) of the metric, 
because this is that component which is necessary and sufficient in order to obtain the 
multipole moment itself, as opposed to some time derivative of it as would be deduced from 
the Oi and ij components (this is very important because we are looking for a term in 
which is a constant). At the linearized approximation, the "far-field" quadrupole moment 
as seen from X + simply reduces to the canonical moment M^- , and from Eq. (j3.16aj) we get 
(with n lj = n l n j — 



G/Ca = 4^ & M l0 {u) + ■■■. (3.20) 



6G 



Here we focus on the term of the form n^r~ 3 f(u), and the ellipsis refer to all the other 
terms, either involving some other multipolarities with i 7^ 2, or a power of 1/r different 
from 3. Consider now the corrections brought about by the non-linear terms to be added to 
the linearized expression ()3.20|) . and write 
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M™ d (u) + ■ ■ ■ , (3.21) 



where the far-field or "radiative" quadrupole moment (i.e., as seen from X + ) is denoted 
by M?-. The non-linear terms in (j3.21|) introduce many couplings between the different 
multipole moments, and there are a lot of possibilities in the general case. However, things 
are much simpler in the case of the BSS owing to the particular structure of the moments 
as determined in Eq. (|3.19|) . Notably, the £-th time derivative of the BSS moment is 
always a constant. Then we find that at 3PN order the most general form of the allowed 
non-linear terms in M-l ad reads as 

CM C 
Mg* = M tJ + 7 -j-Mij + e - M {i M 3) 

G G G 

+ a- M k{l M j)k + 4>- M k M ijk + 6 — M k M ijk 

G 2 M 2 G 2 M ( 1 \ 

+ P -^T-*« + V — MM + O J , (3.22) 

where 7, e, a, <f), 6, p, 77 represent some unknown (for the moment at least) coefficients, which 
are in general constant but as we shall see which can also depend on the logarithm of the 
distance r. Note that Eq. ()3.22)l might a priori contain also some non-local contributions, say 
/_ dv MMij (v) , or (worse) dv M 2 Mij(v), but we shall discuss such terms in Appendix 
lAl and show that they do not appear. The most important terms for the present purpose 
are the two last ones, with coefficients p and 77, which involve the cubic-order multipole 
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interactions M 2 x My and M x Mj x Mf (where M is the mass and Mi the mass dipole of 
theBSS). 

The coefficient we are looking for in the 3PN BSS source-type quadrupole moment iy 
was denoted by C, and we now write the corresponding term as 

8 c I ij = C^-j-V< i V». (3.23) 

The other terms in the BSS quadrupole have a different structure which has already been 
displayed in Eq. (|3.1 2|) . The result C = 4/7 of Section IIII Al will be recovered by the 
following method. As we noticed after Eq. (|3.18jh the coefficient C is necessarily the same 
for the source-type and canonical-type moments, hence 

c 2 m 3 

5cM i3 = C V {i V j) . (3.24) 

Consider next the radiative-type moment (|3.22|) . and look for the modification of C induced 
by non-linearities. In order to do this we have to remember the fact (see Ref. j38|) that 
both the canonical and source moments Ml and II for general matter systems admit a PN 
expansion which is "even" up to 2PN order, with the first "odd" correction being at the 
2.5PN level. This simple fact immediately shows that it is impossible that the "odd" terms 
shown in Eq. (|3.22)1 . which carry explicitly in front the odd powers 1/c 3 and 1/c 5 , contribute 
to a term at the 3PN order. So we conclude that the sought-for modification of C can come 
only from the two last terms, with coefficients p and rj. Taking into account the Newtonian 
results My = 2MV< i V j > + 0{c~ 2 ) and M i = MV i + 0{c' 2 ) we find 

C 2 M 3 

b c Mlf = (C + 2p + V ) . (3.25) 



(Here our notation <5 C M™ means that we are considering the complete term having the 
above indicated structure.) 

Next we come to the central part of this investigation, namely the computation of the 
two coefficients p and rj. This task is not so easy because p and i] are in factor of some 
cubically non-linear terms. We have obtained them by straightforward application of the 
MPM algorithm of Refs. 0, 33[. Actually the value of p, corresponding to the interaction 



M 2 x Mij, is already contained in the result of the calculation of the gravitational wave 
"tails-of-tails" in Ref. [?} . For convenience we relegate the details of this non- linear iteration 
to Appendix |XJ and simply quote here our end results: 
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1271 
9 ~ 735" 

2918 
71 = 735" 

For completeness we give also the known values of three other constants in 1)3.22)1 : 

7 7 20 . 

7 = -, €=-, a = —. 3.27 
' 2 ' 3 ' 21 v ; 

These values come from Eq . (2.8a) in Ref. 0] for 7, Eq. ()A9aj) in Appendix A below for e, 
and the Table 2 in Ref. |33j | for a . (We have not computed the coefficients <fi and 9.) 
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As we see from (J3.26|) . both p and r] depend on the logarithm of r/ro, where tq is the same 
constant as in the source multipole moments (j2.14j) . but we nicely find that the logarithms 
cancel out in the relevant combination of these coefficients which enters Eq. (j3.25|) . In 
fact, we can argue that the cancellation of the logarithms must necessarily occur because 
nowhere in the far-zone expansion of the BSS metric (|3.8|) - (|3.9|) can such logarithms of r/r 
be generated. Hence we get 



Let us now check that there are no gauge effects, linked to the non-geometrical nature 
of our definitions for the multipole moments, concerning the particular term we consider 
in (J3.28|) . in the sense that the coordinate transformation between the canonical metric 
coefficient h^ n and the corresponding BSS one, computed from the expansion of Eqs. (13. 8 j) - 
()3.9|) and given by A4(h 00 ) in the MPM formalism, has no effect on this particular term. The 
proof goes by noticing first that the coordinate transformation at the linearized level is given 
by the gauge vector (j2.9|) parametrized by the four source-type moments Wl, Xl, Yl and 
Zl- The latter moments were given in Eqs. (5.17)-(5.20) of Ref. 0, and we have provided 
in Eqs. ()2.30|) above their new forms in terms of surface integrals. An important point is 
that the moments Wl, X l , Y l , Z l have been defined in such a way that they admit some 
non-zero finite limits when c — > +00; in other words, they "start at Newtonian order" and 
their Newtonian limit is non-zero. By using this fact together with dimensional analysis, it 
is a simple matter to write down their structures in the case of the BSS, in a manner similar 
to what we did for the moment Ml in Eq. ()3.19|) . We find 



while Zl is a current-type moment so Zl = with our choice of origin for the BSS. The 
structures ()3.29|) imply that the potentially dangerous term, which is proportional to M 3 , 
must necessarily appear in these moments at order 1/c 6 relatively to the Newtonian order. 
Next one readily shows, again on dimensional grounds, that the only possible modifications 
of the quadrupole moment My in (|3.2(J|) which are due to the gauge transformation, take the 
forms Wij/c 2 , Xij/c 4 or Y^jc 2 . It is therefore impossible, because of the latter extra factors 
1/c 2 or 1/c 4 , that a dangerous term be generated in this way at the 3PN order. Similar 
arguments are even more easily applied at non-linear order in the coordinate transformation: 
For instance we find that it is impossible that a non-linear coupling of the type Ml x Wp, 
or Wl x X P , has the correct structure at 3PN order in the quadrupole moment. 

This check being done, we conclude that the coefficient we predicted for the relevant 
term in Eq. ()3.28j) represents exactly the quadrupolar contribution in the retarded far-zone 
expansion (at X + ) of the BSS metric (and by the reasoning of Section HII Al we know that we 
can use the BSS metric in standard harmonic coordinates). Now the point is that we can 
also compute directly this contribution in the far-zone expansion of the BSS metric by using 
the formulas ()3.8|) - (j3.9|) . Since the far-zone expansion is to be done at retarded time (and 
not at time t = const as we did in Section Ull Aj) . we must for this calculation substitute t 
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structure as the one of M L given by (|3.19|) | , (3.29a) 
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by u + r/c in Eqs. (13.10(1 and only afterwards take the limit r — > +00 (holding u = const). 
In this way we obtain the quadrupolar piece oc n y '/r 3 in the 00 component of the metric, 
and as we have proved before we are allowed to identify the term therein having the correct 
structure with the one computed in (|3.28() . A simple Mathematica calculation reveals that 
the term in question has the coefficient: 8. Therefore C + y = 8 and we obtain 

C = ^ , (3.30) 

in complete agreement with our previous finding ()3.12j) . and in support of the value for the 
kinetic ambiguity parameter: ( = —7/33. 



APPENDIX A: NON-LINEAR MULTIPOLE INTERACTIONS 

This Appendix is devoted to the computation of the cubically non-linear coefficients p and 
7] entering Eq. ()3.22j) . Following the MPM algorithm of Ref. |6J we must first compute the 
non-linear cubic source term, say A^ n3 , which is composed of the sum of a quadratic-order 
piece made out of products between h can i and h can 2 (and of course their gradients), and 
of a purely cubic-order piece, involving three factors h ca , n \. The linearized metric has been 
given in (|3.16j) . and from it all subsequent iterations are generated by the MPM algorithm, 
which is the same as in Eqs. (I2.11|) - (j2.12j) except that, for the "canonical" construction, the 
gauge vector is <f1 = and we use M L and Sl as moments instead of II, J^. The cubic 
source term A^ n3 [/i can i, ^can2] is inverted by means of the retarded d'Alembertian operator, 
regularized by the specific finite part FPb=o of Eq. ()2.12j) . 
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The metric at cubic order reads then 

^can3 = M can3 + "^canS > (A2) 

where ^^3 represents a particular homogeneous solution of the wave equation, such that 
the harmonicity condition duh^^ = is satisfied. Below we shall not need to compute 
t>can3; since we shall simply have to invoke the fact that its 00 component, v^ n3 , is made of 
multipolarities I = or 1 only (see e.g. Eq. (2.12a) in |33j), and will thus always be zero for 
the quadrupole case £ = 2 of concern to us here. 

Consider first the cubic interaction M x M x My. In this case A^ n3 is already known 
from previous work, Eq. (4.16a) of Ref. 0, which gives, for the needed 00 component, 

A°° n3 = n ab M 2 [-516 r' 7 M ab {u) - 516r- 6 M a6 ( M ) - 304 r" 5 M ab (u)j . (A3) 



To integrate we use the formulas given in the Appendix A of Ref. |33[ . We limit ourselves 
to the computation of the part u^ n3 of the algorithm since v^ n3 = for the particular 
multipole interaction we consider. In all this calculation we use the fact that the second 



7 In this Appendix we pose G 
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time derivative of the BSS quadrupole moment is constant, hence M%j = 0. First of all, by 
straightforward use of the integration formula (A. 16a) of [33] we obtain 
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The last term is a priori more delicate in the case of the BSS because we know from 
Eq. (A. 13) of [3^ that it could generate an integral depending on the whole time evolution 
of the system, such as a "tail" integral of the type J{u) = f_ dv M 2 Mij(v). The value 
of an integral like J would be quite ambiguous within our MPM-based approach. Indeed, 
on the one hand, the boosted system that we finally consider has My-(u) proportional to 
u 2 , so that Mij(u) vanishes identically. We would therefore expect, from this point of view, 
that J, being the integral of a vanishing integrand, is zero: J = 0. On the other hand, we 
can perform the w-integral in J to get J = M 2 [My(u) — Mjj-(— 00)]. Now, as we recalled 
above, the MPM framework on which we base our discussion assumes that we start initially 
by considering systems such that the multipole moments become time-independent in the 
remote past (before some finite instant —T). For such systems, the second contribution 
in the latter expression for J vanishes, and we would get: J = M 2 Mij(u), which does not 
vanish in the case of a boosted source. 

However, this ambiguous situation does not appear in the present calculation. Indeed, 
the logic of our calculation is the following. To derive the MPM result flA4|) we had to 
initially assume that My- (14) tends fast enough toward a constant in the infinite past. Then, 
after having done the MPM iteration we get the form (jA4j) . Now, starting from the ex- 
plicit expression (jA4|) we want to relax the original assumption about My(u), and consider 
a deformation process in which M^-(ti) interpolates between an initial MPM-like Mij(u) 
(tending to a constant in the past) and a final (u) of the form of u 2 . The question 
is then to know (i) whether the R.H.S. of Eq. (|A4|) admits any limit after this continuous 
deformation process, and (ii) what is the value of this limit. In mathematical terms the 
question is essentially a question of interchange of a limit with an integral operation, i.e. 
whether lim n f f n (u)du = f\im n f n (u) du holds, for a sequence of functions f n (u). As we 
know the answer is positive under "good" conditions, for instance of uniform convergence, 
or more generally (Lebesgue's theorem) of dominated convergence, which says essentially 
that if \ f n (u)\ < g(u) and if J g{u)du is finite, then we can interchange the limits. In our 
case, we can use Lebesgue's theorem of dominated convergence (see, e.g., the book [39]), 
which is both simple and powerful. 

The only delicate term in the R.H.S. of (jA4|) is the last, integral term. We must interpolate 
between some initial My(w) which vanishes in the past to become a constant in the future, 
and a final Mfj (u) which is always constant. It is clear that we can do this interpolation in 
a way that |My(«)| remains always bounded. Using such a bound in the last term of (jA4|) . 
which is a three-dimensional (retarded) integral, we easily see that, under the assumption 
of a bounded |My(«)|, the integrand is bounded by a (positive) function which is mtegrable 
(because of the fast convergence brought by the r~ 5 factor, together with the 1/r factor 
contained in the propagator). Therefore, we can indeed interchange the limiting process 
and the integration one, and conclude that the limit of the L.H.S. of Eq. (|A4|) exists, and is 
simply given by replacing Mij(u) in the R.H.S. by its limiting expression for a BSS which 
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is: M^ ss (w) proportional to u 2 . As, under this limit, the integrand of the last term in (|A4|) 
becomes time- independent, we can explicitly compute the limit by replacing the retarded 
propagator by a Poisson integral. Hence, we get 
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One easily checks that the result ()A5|) agrees with the more formal way of doing the cal- 
culation which consists of applying the formula (A. 13) of Ref. |33| to the BSS case. As we 
can see, the last term in the R.H.S. of the latter formula, which is given by a "tail" integral, 
vanishes when we insert the BSS quadrupole into the source term of the retarded integral, 
since this source term involves in fact the fifth time derivative of the BSS quad rupole which 
is zero. On the other hand, the logarithmic term in the formula (A. 13) of |33j does remain, 
and then we recover exactly Eq. (|A5|) — with the same constant tq on both sides of the 
equation. The proof that we have detailed above, based on Lebesgue's theorem of domi- 
nated convergence, rigorously justifies (for the case at hand) that one can "blindly" use the 
formulas in the Appendix of [33j even for the case of the BSS. 
Finally, gathering Eqs. (|A4|) and (|A5|) we obtain 
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which shows by comparing to (j3.21|) - ()3.22|) that the coefficient we are seeking is 
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The cubic multipole interaction M x Mi x Mj is much longer to obtain because we are 
obliged to compute it from scratch (no earlier results in the literature are available). At 
linearized order the metric reads 
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Straightforward calculations following the MPM algorithm then yield the part of the metric 
at quadratic order corresponding to the multipole couplings MxMj and Mj x Mj, 
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(A9b) 
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Using such expressions (1A8J) and (|A9Jl we next obtain the source term at the cubic-order 
approximation M x Mi x Mj. We are interested only in its 00 component which is then 
found to be 
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The integration proceeds exactly in the same way as in Eqs. ()A4|) - (jA5|) . Again, the problem 
is that of interchanging a limiting process Mi (it) — > Mf (u) with the retarded integration, 
and this can be proved by using the Lebesgue theorem, because the various powers of 1/r in 
the R.H.S. of (jA10|) ensure that, in a process where |Mj(m)| remains bounded, the integral 
is bounded by a positive convergent integral. And again the result ag rees with the one we 
would formally obtain by using the formulas in the Appendix of |33j]. As before we have 
t>can3 = from the argument concerning the multipolarity i = 0, 1 of this special piece. Our 
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result is then 
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from which one recognizes on comparison with (j3.21j) - ()3.22j) that 
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